Spectral properties and pattern selection in fractal growth networks 
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A model for the generation of fractal growth networks in Euclidean spaces of arbitrary dimension 
is presented. These networks are considered as the spatial support of reaction-diffusion and pattern 
formation processes. The local dynamics at the nodes of a fractal growth network is given by 
a nonlinear map, giving raise to a coupled map system. The coupling is described by a matrix 
whose eigenvectors constitute a basis on which spatial patterns on fractal growth networks can 
be expressed by linear combination. The spectrum of eigenvalues the coupling matrix exhibits a 
nonuniform distribution that is reflected in the presence of gaps or niches in the boundaries of 
stability of the synchronized states on the space of parameters of the system. These gaps allow for 
the selection of specific spatial patterns by appropriately varying the parameters of the system. 
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I. INTRODUCTION 

Pattern formation processes in regimes far from equi- 
librium often take place on media that are nonuniform 
at some length scales. The nonuniformity may be due 
to the intrinsic heterogeneous nature of the substratum, 
typical of pattern formation in biological contexts, or 
it may arise from random imperfections or fluctuations 
in the medium. Such heterogeneities can have signifi- 
cant effects on the form of spatial patterns, for exam- 
ple, they can induce reberverators in excitable media and 
defects can serve as nucleation sites for domain growth 
processes. Recently, there has been much interest in the 
study of spatiotemporal dynamical processes on nonuni- 
form or complex networks. In this context, coupled map 
lattices 0] have provided fruitful and computationally 
efficient models for the investigation of a variety of dy- 
namical processes in spatially distributed systems. In 
particular, the discrete-space character of coupled map 
systems makes them specially appropriate for the inves- 
tigation of spatiotemporal dynamics on nonuniform net- 
works that can represent models of heterogeneous me- 
dia. Phenomena such as pattern formation, spatiotem- 
poral intcrmittcncy, nontrivial collective behavior, syn- 
chronization, phase-ordering, etc., have been extensively 
studied in coupled map systems defined on fractal lattices 
13,01, hierarchical structures p|, trees 0, random graphs 
[g , small- world networks , and scale-free networks @ . 

An especially interesting class of nonuniform geome- 
tries comprises fractal growth networks whose branch- 
ing structure, self-similar scaling features and lack of 
translation symmetry can give rise to several distinct 
characteristics in both their dynamical and spatial prop- 
erties. Fractal growth structures appear in nonequilib- 
rium growth processes which are common in many areas. 
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Examples of such phenomena include diffusion-limited 
aggregation, dendritic solidification, bacterial growth, 
viscous fingering, capillarity, electrodeposition, Lapla- 
cian growth problems, etc. flO| . 

In this article we consider discrete reaction-diffusion 
processes and pattern formation taking place on fractal 
growth networks. The spatiotemporal dynamics corre- 
sponds to a coupled map system defined on the geometri- 
cal support of a fractal growth network. Although many 
growth structures found in nature have random features, 
here we study the case of simple, deterministic fractal 
growth networks. We focus on the changes occurred in 
spatiotemporal patterns as a result of the fractal-growth 
connectivity that describes the interactions in the sys- 
tem. In Sec. II, a model for the construction of fractal 
growth networks in any Euclidean space and a general no- 
tation for their treatment are introduced. In Sec. Ill, the 
coupled map models defined on fractal growth networks 
are presented. The diffusion coupling among neighbor- 
ing sites on a network is described by a matrix. The 
spectrum of eigenvalues and eigenvectors of the coupling 
matrix is analyzed in Sec. IV. The eigenvectors consti- 
tute a complete basis on which spatial patterns can be 
expressed by a linear combination. Section V contains a 
study of the stability of spatially uniform, periodic pat- 
terns on a fractal growth network for a local dynamics 
given by the logistic map. Distinct features, emerging 
as a consequence of the fractal nature of the networks, 
allow for the selection of specific spatial patterns as the 
parameters of the system are changed. Conclusions are 
presented in Sec. V. 



II. MODEL FOR FRACTAL GROWTH 
NETWORKS 

Fractal growth networks can be generated in any Eu- 
clidean space of dimension d by generalizing the frac- 
tal growth model of Vicsek [||. Starting from a d- 
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dimensional hypercube, it is divided into 3 equal hy- 
percubes. Of these, only the central hypercube plus the 
2 d hypercubes connected to its vertices are kept. This 
process is repeated for each of the 2 d + 1 resulting hy- 
percubes. At a level of construction L, the network con- 
sists of N = (2 d + 1) L hypercubical cells or nodes whose 
coordinates can be specified by a sequence of symbols 
(a\a.2 ■ ■ ■ where ai can take any value in a collection 
of 2 d + 1 different digits forming an enumeration system 
in base (2 d + 1), which we denote by {ei, e 2 . . . , £2 d + 1 +i} • 
For example, for d = 2, we may chose at £ {0, 1, 2, 3, 4}. 

To fix the labels, we start at level L = 1 where the 
labels of the 2 d + 1 cells in the network have the form 
(ail). The label (0) is assigned to the central cell, and 
the others 2 d cells are labeled by (ai) ^ (0) in such a 
way that the labels (o^) and (a>j) of the two cells located 
at the opposite ends of each diagonal passing through 
the central cell satisfy mod 2 d +1 (ai + (Xj) — 0, where the 
addition at + otj is defined modulo (2 d + 1). Figures 1(a) 
and 1(b) show this notation for the fractal growth net- 
works embedded in Euclidean spaces of dimension d = 2 
and d = 3, respectively. 
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FIG. 1: (a) Fractal growth network embedded in Euclidean 
space of dimension d = 2 at levels of construction L = 1, 
L = 2 and L = 3, showing labels on the cells, (b) Fractal 
growth network embedded in dimension d = 3 at level of 
construction L — 2, showing labels on the cells. 



When the fractal network grows from level of construc- 
tion L to level (L + 1), each cell {a\a 2 ■ ■ ■ chl) is subdi- 
vided into 2 d + 1 cells, scaled down by a longitudinal fac- 
tor of 3, and which are now labeled by (a±a2 . . . a^ai+i), 
where the first L symbols of the sequence, a±a2 . . .cv_l, 
are the same as in the parent cell, and to distinguish the 
(2 d + 1) daughter cells, the last symbol is assigned 

similarly to the labels (a±) at the construction level 
L = 1. From the construction process, it follows that the 
fractal dimension of the network is dt = ln(2 d + l)/ln3. 

The set of cells connected to the cell labeled by 
(cvia2 • • -(Xl) at a level of construction L is defined as 
the neighborhood (nearest neighbors) of this cell, and 
will be denoted by J\[(a\a 2 . . . ah)- Three types of cells 
can be identified in a fractal growth network embed- 
ded in a d-dimensional Euclidean space: (i) centers, that 
are connected to 2 d + 1 other cells; (ii) joints, that are 
connected to 2 cells; and (iii) edges, connected to just 
one cell. A sequence (aio^ . . . oll) can be written as 
(pt\ . . . otL-sOi s L _ s+ i) for some s £ {1, 2, . . . ,L — 1}, and 
where af means the sequence of s symbols ev^. A cen- 
ter cell is labeled by a sequence (evicts . . . aL- s a s £~_ L s+1 0). 
Its neighborhood set consists of 2 d cells, labeled by the 
sequences (aia 2 ■ ■ ■ ajl 1 s+1 l), (aia 2 ■ ■ ■ ai_ s a^l 1 s+1 2), 
. . ., and (ctict2 ■ ■ ■ oiL^ s a s £[} a+l 2 d ) , respectively. A joint 
cell labeled by (a± . . . aL~sCt s L _ s+1 ), with a^-s+i 7^ 0, 
has a neighborhood set with 2 elements. One is the center 
cell labeled by (a\ . . . oi^_ s a s L [} s+l Q). The other neighbor 
is another joint cell labeled by (a\ . . . aL- s -ictL-s+i/3 s ), 
if cth-s = 0; or by («i . . . a L - s -iO(3 s ), if a L - s ^ 0; where 
P is an allowed symbol that satisfies mod 2 ii + i(aL- s +i + 
0) = 0. This procedure identifies any two joint cells as 
reciprocal neighbors. On the other hand, an edge cell 
labeled by the sequence (a\ . . . cvL_ s a^_ s+1 ) can not be 
assigned a reciprocal neighbor from the above procedure; 
thus this edge possesses only one neighbor, the center cell 
labeled by («i . . . a i _ s a£l 1 s+1 0). 

Because the symbols ai belong to an enumeration sys- 
tem in base (2 d + 1), a cell in the network labeled by 
the sequence (a±a2 ■■ - oil) can be univocally associated 
to an integer index i = 0, 1, 2, . . . , (2 d + 1) L — 1, by the 
rule 

L 

(a 1 a2...a L )^i = Y J ^{2 d + \) L - ] . (1) 

7=1 

As an illustration, consider the fractal growth network 
embedded in an Euclidean space of dimension d — 2, as 
in Fig. 1(a). At the level of construction L = 3, the cell 
labeled by (233) = (23 2 ) is a joint cell and has index 
i = 68, according to the rule Eq. QJ. Its two neighbors, 
labeled by (230) (a center cell) and (022 = (02 2 ) (an- 
other joint cell), are assigned indexes i — 65 and i = 12, 
respectively. Similarly, the cell (422) = (42 2 ) is an edge 
with associated index i — 112; its only neighbor is the 
center cell (420), corresponding to i — 110. 
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III. COUPLED MAP LATTICE MODEL FOR 
PATTERN FORMATION 

A network might be considered as the spatial support 
of a dynamical spatiotemporal processes with either dis- 
crete or continuous time. Here we consider reaction- 
diffusion and pattern formation phenomena on fractal 
growth networks. By associating a nonlinear function 
to each cell of a given fractal network and coupling these 
functions through nearest-neighbor diffusion interaction, 
we define a coupled map lattice that describes a reaction- 
diffusion dynamics as follows 

z t+1 (i) = /(z t (i))+ 7 J2 fatf) - **(*)) . ( 2 ) 

where xt(i) represents the state of the cell having index 
i, assigned by Eq. at discrete time t; f(x t (i)) is a 
nonlinear function specifying the local dynamics; is 
the neighborhood set of the cell with index i; and 7 is a 
parameter expressing the coupling strength among neigh- 
boring cells and it plays the role of a homogeneous diffu- 
sion constant. The form of the coupling term in Eq. Q 
is usually called backward diffusive coupling and corre- 
sponds to a discrete version of the Laplacian in reaction- 
diffusion equations. This coupled map model can be gen- 
eralized to include other coupling schemes, nonuniform 
coupling, or continuous-time local dynamics. 

Equations (0) can be written in a vector form as 

x i+ i = f(x t )+ 7 Mx t , (3) 

The state vector x t has N components [x t ]i = x t (i), 
i = 0, ...,N— 1, corresponding to the states Xt{i) — 
Xt{a.\OL2 ■ ■ ■ ol£) of the cells on the network. The N x N 
real, symmetric matrix M expresses the coupling among 
the N components Xt(i). For a fractal growth net- 
work embedded in an Euclidean space of dimension d 
at the level of construction L, the components of the 
corresponding matrix M, denoted by M(i,j) = 
0,1,..., N - 1), are 

fl, if jeJV(i) 

M(iJ) = M(j, i) = l- , if i = j (4) 

1 0, elsewhere, 

where is the cardinality of the neighborhood set 

Af(i) — Af(aiOt2 ■ ■ - oil)- The matrix M plays the role 
of the spatially discrete diffusion operator on these net- 
works, similar to the Laplacian in a spatially continuous 
reaction-diffusion equation. 

IV. SPECTRUM OF THE COUPLING MATRIX 

The spatial patterns that can take place on fractal 
growth networks are determined by the eigenmodes of 
the coupling matrix M, similarly to reaction-diffusion 
processes on regular Euclidean lattices [Til . On the 



other hand, the stability of the synchronized states is 
related to the set of eigenvalues of M. 

In order to analyze the eigenvector problem, consider a 
fractal growth network embedded in an Euclidean space 
of dimension d, at level of construction L, on which a 
spatiotemporal dynamics has been defined in the vector 
form of Eq. © . The complete set of orthonormal eigen- 
vectors of the corresponding matrix M can be described 
as the superposition of two distinct subsets of eigen- 
modes. One subset, which will be denoted by {u m „}, 
contains those eigenvectors associated to non-degenerate 
eigenvalues; and the other subset comprises the eigenvec- 
tors corresponding to degenerate eigenvalues of M, and 
will be represented by {v^ n } (the indices refer to the de- 
generacy, as explained bellow). Thus, the complete set 
of eigenvectors of M is {u mn } U {vj^}. Each eigenvec- 
tor describes a basic spatial pattern that may arise on 
a fractal growth network characterized by an embedding 
dimension d and a level of construction L. As illustra- 
tion, we shall show the eigenvalues and eigenvectors of 
the coupling matrix corresponding to the fractal growth 
network embedded in the plane. 

Non-degenerate eigenvectors 

The eigenvectors of M belonging to the non-degenerate 
subset {u mn } satisfy 

ivr — h J to = 0; n = 

!VlU m „ - mn U mn , I m = 1 )i; n = lf j gm-l. 

(5) 

where b mn is a eigenvalue associated to the eigenvec- 
tor u mn . The i component of vector u mn is [u mn ]j = 
Umn(aiC2 • ■ ■ according to the rule in Eq. QJ. The 
number of eigenvectors in this subset and that of their 
corresponding eigenvalues, denoted by v(b), is 

u(b)=l+J2^- 1 = ^ 1 - (6) 

m—l 

Any eigenvector u mn in the non-degenerate subset {u mn } 
is characterized by the following property: its compo- 
nents corresponding to cells of the same type on the 
network and located at the same distance (measured in 
number of cells) from the center cell (0 L ) are identical. 
Because of this property, we also refer to the elements in 
the subset {u mn } as symmetric eigenvectors. 

The index m corresponds to the level of construction at 
which the eigenvalue b mn first appears, and it is related to 
the length of the longest diagonals that arise as the frac- 
tal growths up to the level L. At a level of construction 
L > 1 there is a total of (2 d + 1) L_1 + 1 diagonals on the 
network, which can be classified in L distinct families or 
groups according to their lengths, measured in terms of 
number of cells on the diagonal. Each different family of 
diagonals can be identified by the index m = 1, 2, . . . ,L, 
and each family contains all the diagonals of the network 
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that have 3 m cells, and thus can support a wavelength 
of 3 m . The index m gives the step in the construction 
process at which the family of diagonals of length 3 m 
has first appeared. For a network defined at level L, 
there are (2 d x (2 d + i^-™- 1 ) diagonals in each fam- 
ily characterized by m = 1, 2 . . . , L — 1, and 2 diagonals 
in the family corresponding to m = L. For example, in 
Fig. 1(a) with d = 2 and L = 3, there are 3 families of 
diagonals according to their lengths: the family m = 1 
contains 20 short diagonals with a length of 3 cells; the 
family m = 2 has 4 medium diagonals measuring 9 cells; 
and the family m = 3 possesses 2 long diagonals with 27 
cells. Additionally, we must count the family m = as- 
sociated to diagonals having 1 element; this corresponds 
to the spatially homogeneous eigenvector. Thus, the in- 
dex m = 0, 1, . . . , L, indicates the level of construction at 
which the family of diagonals characterized by the index 
m and having the same length of 3 m cells has appeared. 
The family of diagonals identified by m remains in the 
fractal as the network grows up to level L. 

The index n counts the number of distinct symmetric 
eigenvectors that have originated each time a new family 
of diagonals of length 3 m cells appear, and it depends 
on the number of intersections that a new diagonal of 
length 3 m has with the other diagonals already present 
in the network. There is one intersection for each 3 cells 
in a diagonal; thus the number of intersections is 3 m_1 . 
These intersections determine the different wavelengths 
that can be formed on a diagonal of length 3 m . The num- 
ber of these different wavelengths having reflection sym- 
metry about the long diagonal of length 3 m originated at 
step m can be counted by n = 1, 2, . . . , 3 m_1 . 

The spatially homogeneous eigenvector of the matrix 
M, which we denote by u o, belongs to the subset {u mn } 
and its components N are 



u o(ai 



= AT-i/2. 



VL, Vat 



(7) 



fractal growth network embedded in an Euclidean space 
of dimension d = 2, at level of construction L = 2. 




FIG. 2: Fractal growth network corresponding to embedding 
dimension d — 2 and level of construction L — 2. (a) The 
five non-degenerated, symmetric eigenvectors u mn and their 
corresponding eigenvalues, (b) Degenerate eigenvectors v mn 
and their corresponding eigenvalues. 

The set of eigenvalues arising from Eq. (J5J may be or- 
dered by decreasing value. By Gershgorin's theorem [l3| . 
the homogeneous eigenvector possesses the largest eigen- 
value of M, which is &oo = 0. In the case of embedding 
dimension d — 2, the smallest eigenvalue for large L is 
found to be 



lim 31.-1 = — (4 - 

L — >oo 



(10) 



and 



For each index m > 1, the eigenvalues 

{^mb b m 2i • • • 5 ^m3 m 

-1} arise in groups of three, b 

run- 

b mn > y b mn n , where 



u 00 = -==001(1, 1, 1) 
V N 



(8) 



Since the eigenvectors of M are mutually orthogonal, all 
others eigenmodes in either subset {u mn } or {v^„} must 
satisfy 



u„ m (ai . . .a L ) = 0; Va fc , m^O; 



E v L(ai--«i) = 0; Va k 



(9) 



that is, the sum over the components of any eigenvector 
of M, different of Uoo, is zero. 

Figure 2(a) shows the v(b) — 5 non-degenerate eigen- 
vectors and their associated eigenvalues of the coupling 
matrix corresponding to a fractal growth network embed- 
ded in dimension d = 2 and having a level of construc- 
tion L = 2. Figure 2(b) shows the subset of degenerate 
eigenvectors v^ s and the eigenvalues corresponding to a 



n' = 2 x 3" 1 " 1 - n + 1; n" = 2 x 3™- 1 + n. (11) 

The sum of the eigenvalues in each of these groups is 
constant, and gives 



-(2 d + 4). 



(12) 



Because of (I12fl . the eigenvalues associated to the non- 
degenerate eigenvectors satisfy 



L 3"- 1 

bmn 

m— n— 



-(2 d + 4) 

m— 1 



= -(2 d - 1 + 2)(3 i -l). (13) 



Figure 3 shows the spectrum of eigenvalues {b mn }, in- 
dicated by empty symbols, for a fractal growth network 
embedded in an Euclidean space of dimension d = 2, at 
successive construction levels L. Eigenvalues associated 
to degenerate eigenvectors of the coupling matrix M, to 
be discussed next, are also shown in Fig. 3. 
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FIG. 3: Spectrum of eigenvalues of the coupling matrix at 
increasing construction levels L, for a fractal growth network 
embedded in dimension d = 2. (a) Eigenvalues b mn are indi- 
cated by empty symbols: & o (□); &n (0); b 2 „ (A); b 3 „ (V); 
bin (0)- Other symbols indicate eigenvalues a mn as follows: 
ain (■); a2n (•); a- in (a); a4n (x); a 5n (+). (b) Magnification 
of the dotted region in (a). 



Degenerated eigenvectors 

The subset of degenerate eigenvectors {v^ (n } of the 
matrix M satisfy 

m=\,2,...,L; n = 1, 2, . . . , 3™- 1 ; 1 J 

where a mn is the eigenvalue associated to a group of f2 
degenerate eigenvectors {v} nn , v r 2 n „, . . . , v^„} belonging 
to the subset {v^ ln }. The index g runs from 1 to a value 
f2 and counts the different eigenvectors associated to the 
degenerate eigenvalue a mn , as it will be shown bellow. 
The integer indices m and n label different eigenvalues 

The i— component of a vector v^ n corresponds to a cell 
of the network labeled by the rule Eq. QJ, i.e., [vf nn ]j = 
v mn( a i Q! 2 ■ ■ • oil)- The eigenmodes in the subset {vf ml } 



are characterized by the following two properties, 

v^„(a 1 ...a L _ m m ) = 0; m = 1, 2, . . . , L - 1; (15) 

that is, all the components of v^ m corresponding to cen- 
ter cells formed in the first L — m levels of construction 
vanish; and 

v a mn (a-La2---OL L -iti) = 0; (16) 

that is, the sum of the components associated to the cen- 
ter cells in a network spatially described by a vector v^ s , 
is zero. 

In contrast to the eigenvectors in the non-degenerate 
subset, an eigenvector v^ s is non-symmetric; however it 
exhibits a partial regularity, having all its components 
corresponding to center cells equal to zero if they origi- 
nated at the level of construction L — m. The index m 
counts the number of additional levels of construction in 
which non-vanishing center cells have appear up to level 
L; and its possible values are m = 1,2, ... ,L. There 
are (2 d + l) L ~ m vanishing center cells in an eigenvec- 
tor v^ s . Each null center cell in the eigenvector is the 
center of 2 d_1 — 1 diagonals, except the cell (0 L ) that 
is the center of one additional diagonal. Thus, adding 
the additional diagonal of the center cell (0 L ), there are 
(2 d + l) i_m (2 d_1 — 1) + 1 diagonals whose centers are 
null center cells. The index n counts the number of pos- 
sible intersections that these diagonals have with other 
diagonals; its possible values are n— 1,2,..., 3 m_1 . 

Because an eigenvector v^ s is a non-symmetric mode, 
the 2 edge cells at the ends of each of the (2 d + 
l) L ~ m (2 d ~ 1 — 1) + 1 diagonals having null center cells 
are independent of each other. Therefore, there exist 
2[(2 d + l) L ~ m {2 d ~ 1 - 1) + 1] edge cells with this charac- 
teristics, for given values of m and n. However, due to the 
orthogonality property, Eq. @ , the total number of lin- 
early independent eigenvectors v^ s possessing the same 
indices m and n is fi = 2[(2 d + l) L - m (2 d - 1 - 1) + 1] - 1. 
The index g counts the number of independent eigenvec- 
tors associated to an eigenvalue a mn , and it may take 
the values g = 1, 2, . . . , CI. In this fashion, the subset of 
degenerate eigenvectors {v^} of the coupling matrix M 
corresponding to a fractal growth network embedded in 
dimension d and at level of construction L is fully de- 
scribed. 

In the case of embedding dimension d = 2, there are 
5 L ~ m vanishing center cells in an eigenvector v^ s and 
there are 5 L_m + l diagonals whose centers are null center 
cells. The number of linearly independent eigenvectors 
v^ s having the same indices m and n is fl = 2(5 L ~ m + 
1) — 1 = 2 X 5 L ~ m + 1. Thus, the index g may take the 
values g =,2,.. .,2 x 5 L " m + 1. 

The number of different eigenvalues, a mn , that belong 
to the spectrum of matrix M for a fractal growth network 
at construction level L, is 

y[a) = £ 3— 1 = ^1 . (17) 

m— 1 
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At each level m > 1, there appear 3 m-1 new eigenval- 
ues a mn . Similarly to the non-degenerate eigenvalues, 
the a mn can be grouped in sets of three eigenvalues that 
satisfy 



+ a n 



= -(2 d + 4), 



where 



2 x 3 m ~ 1 



1 



2 x 3 r ' 



(18) 



(19) 



and therefore the total sum of the degenerate eigenval- 
ues a mn of a matrix M corresponding to a fractal growth 
network embedded in a dimension d and at level of con- 
struction L gives 



L 3" 



EE 

m—l n=l 



-(2 d + 4) 3m_1 

m—l 

-(2 d - 1 + 2)(3 i - 1 -l). (20) 



Figure 3 shows the spectrum of eigenvalues {a mn } U 
{bmn} for a fractal growth network embedded in Eu- 
clidean space of dimension d = 2 at successive levels 
of construction L. The distribution of the spectrum of 
eigenvalues of the coupling matrix M can be seen as a 
function of the growth process of the fractal structure. 
Note that the full spectrum {a m „}U{6 mn } is always con- 
tained between the eigenvalues 6oo = and b L3 L~i. 

From Eqs. © and Ijl7|l . the total number of distinct 
eigenvalues of M, including both types a mn and b m n, and 
denoted by ^(M) is given by 



v(M) = v{b) + v{a) 



(21) 



Since there appear 3 m_1 new eigenvalues of type a mn at 
each level of construction m, and there are Jl = 2[(2 rf + 
l) i ~ m (2 !i ~ 1 - 1) + 1] — 1 eigenvectors associated to 
each eigenvalue a mn , the total number of independent 

}is 



eigenvectors in the subset {v*^ 



L 

E 

m—l 



n(3 m_1 ) = (2 d + l) J 



1 



(22) 



From Eq. © we know that the number of independent 
eigenvectors in the subset {u m „} is (3 L + 1) /2. Therefore, 
the total number of independent eigenvectors of M is 



1 



1) 



1 



(2 d + 1) L = N, (23) 



as should be expected. 

Figure 4(a) shows the complete spectrum of eigenval- 
ues of M and the degeneracy fraction of each eigen- 
value, for a fractal growth network embedded in an Eu- 
clidean space of dimension d — 2, at level of construction 
L = 5. The degeneracy SI = 2 x 5 5_m + 1 of each of the 
v(a) — (3 5 — l)/2 = 121 eigenvalues a mn is plotted as a 
vertical bar, while the v(b) = (3 5 + l)/2 = 122 different 



eigenvalues b mn are indicated by plus symbols (+) and 
they are non-degenerate. It is clear that both the distri- 
bution of eigenvalues and their degeneracies are nonuni- 
form. The scaling properties of the spectrum of eigen- 
values of the coupling matrix can also be conveniently 
represented by plotting the accumulated sum of the de- 
generacies of all eigenvalues, that is, the measure of the 
spectrum of M (denoted by p) , on the eigenvalue axis for 
large L, as in Fig. 4(b). The resulting graph exhibits the 
features of a devil staircase, a fractal curve arising in a 
variety of nonlinear phenomena. 




l o 



bmn i ^mn 



FIG. 4: (a) Distribution and degeneracies of the spectrum 
of eigenvalues of the coupling matrix M for a fractal growth 
network embedded in d = 2 at level of construction L — 5. 
The eigenvalues 6 mn are indicated with (+) just below the 
zero line, for clarity. The vertical axis shows the degeneracy 
of the eigenvalues {a mn } divided by N, indicated by a vertical 
bar at each eigenvalue, (b) The measure of the full set of 
eigenvalues of M. 

The eigenvectors of the coupling matrix reflect the 
topology of the fractal growth network and they are anal- 
ogous to the Fourier eigenmodes appearing in regular Eu- 
clidean lattices. In this sense, the symmetry properties 
of the non-degenerate eigenvectors {u m „} and the condi- 
tions given in Eqs. Ijl5|l - (|16JI f° r the degenerate eigenvec- 
tors {v^ 5 } represent different wavelengths on a fractal 
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growth network embedded in an Euclidean space of di- 
mension d at level of construction L. 



V. SYNCHRONIZED STATES AND PATTERN 
SELECTION 

Synchronized states in spatiotemporal systems are rel- 
evant since we are often interested in mechanisms by 
which a spatial pattern can be selected in a uniform sys- 
tem that breaks its symmetry as a parameter is changed. 

Consider spatially synchronized, period- K states such 
as Xt{a\...aL) = Xk, V(cti . . . ai); where Xk, (k = 
1,2, ... , K), is a period- if orbit of the the local map, 
satisfying f^ K \xk) = x~k- The linear stability analysis 
of periodic, synchronized states in coupled map lattices 
is carried out by the diagonalization of the matrix M in 
Eq. Js3, and it leads to the conditions 

K 

ni/'(S*)+7Ml<l, ( 24 ) 
fe=i 

where n is any of the i/(M) = 3 L eigenvalues, associ- 
ated to either subset of eigenvectors {u mn } or {v^ s }, of 
the coupling matrix M corresponding to a fractal growth 
network embedded in a d-dimensional Euclidean space at 
level of construction L. 

The nonuniform distribution of the eigenvalue spec- 
trum is manifested in the stability of the synchronized 
states through Eq. (|2"£jl and gives rise to marked dif- 
ferences when compared, for instance, with the bifur- 
cation structure on regular lattices. As an application, 
consider a local dynamics described by the logistic map, 
f(x) = Ax(l — x). In this case, the stability conditions, 
Eq. (|24H . for the period K — 2 P , synchronized state yield 
the set of boundary curves 

S*0*) = n[A(l-2x fc )+7M]=±l- (25) 
fc=i 

For each sign, Eqs. i|25[l give 3 L boundary curves in the 
plane (7, A), corresponding to the different values of fx; 
these curves determine the stability regions of the period- 
2 P , synchronized states on the network. 

The scaling structure for the period-2 p , synchronized 
states in fractal growth networks is similar to that of 
a any lattice described by a diffusive coupling matrix, 
since the form of Eq. (|25[1 is the same in any case. As for 
regular Euclidean lattices 0, > the stability regions 
for the period- 2 P , synchronized states in the (7, A) plane 
scale as A ~ 5~ p , and 7 ~ a~ p , where S — 4.669... 
and a = —2.502 . . . are Feigenbaum's scaling constants 
for the period doubling transition to chaos. However, 
the specific structure of the eigenvalue spectrum of the 
coupling matrix determines the shapes and gaps of the 
regions of stability of synchronized, periodic states. 

From Eqs. 1)25(1 . the boundary curves for the synchro- 
nized, fixed point state (p = 0) on the parameter plane 



(7, A) are given by the straight lines 

A = M7 + 1, A = p 7 + 3; (26) 

which are first crossed for the most negative eigenvalue, 
fi = b L3 L-i. Similarly, the boundaries period-two (p = 
1), synchronized state on a fractal growth network em- 
bedded in an Euclidean space of dimension d = 2 and at 
construction level L = 3 are given by the two sets 

SKomn) = -X 2 + 2\ + 4 + ~/a mn (ja mn -2) = ±1 , (27) 

and 

Sl(b mn ) = -A 2 + 2A + 4 + 7 6 m „( 7 6 m „-2) = ±1;. (28) 

Figure 5(a) shows the boundary curves Eqs. 127ll - 128|) 
on the plane (7, A). The boundary between the synchro- 
nized, fixed point state and the synchronized period-two 
state occurs at A = 3. The upper boundaries (cor- 
responding to —1 in the r.h.s of Eqs. I|27|) - (|28[l ) have 
minima A m i n = 1 + v5 at values 7 m i n = l/b mn and 
7min = l/a mn (for any period 2 P , A min depends on p). 
Figure 5(b) shows a magnification of Fig. 5(a) around 
the minima of the upper boundaries. The distribution of 
the minima 7 m i n and the presence of nonuniformly dis- 
tributed gaps in the boundary curves is a manifestation 
of the nonuniform structure of the eigenvalue spectrum. 
Since the nonuniformity in the distribution of eigenvalues 
persists at any construction level L of a fractal growth 
network, this property allows for regions of stability of 
the synchronized states or gaps characteristic of fractal 
networks and which are not present in other geometries, 
for example in regular lattices, where the distribution of 
eigenvalues of the coupling matrix is uniform and contin- 
uous in the limit of infinite size lattices. 

The set of eigenvectors {u m „} U {vf nn } of the coupling 
matrix M constitute a complete basis (normal modes) 
on which a state x t of the system can be represented as 
a linear combination of these vectors. Thus, the evo- 
lution of x t reflects the stabilities of the normal modes. 
Figure 5(b) shows how a spatially inhomogeneous pat- 
tern may be selected as the synchronized state becomes 
unstable through crossing of the upper boundary; the 
first boundary segment crossed determines the character 
of the instability and the properties of the selected pat- 
tern. For example, consider an initial state consisting of 
a small perturbation of the synchronized, period-2 state 
at parameter values just beyond the boundary segment 
corresponding to the eigenvalue 036, indicated by a cross 
in Fig. 5(b), where this initial state is unstable. The inho- 
mogeneous period-4 final spatial pattern is represented in 
Fig. 6; it corresponds to a linear combination of the three 
eigenmodes {vf 6 ; g — 1,2,3} associated to the degener- 
ate eigenvalue 036 of the matrix M corresponding to the 
fractal network embedded in dimension d — 2 at level of 
construction L = 3. All other modes are unstable in this 
region of parameter space. For any construction level L 
of the fractal network, and any period 2 P , the boundary 
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curve 53(036) = — 1 separates a gap of the synchronized 
state from the stable region for the eigenmodes v| 6 cor- 
responding to 036- Therefore, a transition between these 
two spatial patterns can always be observed in the ap- 
propriate regions of the parameter plane (7, A). 
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FIG. 5: The boundary curves S3 = ±1 given by Eqs. 
127I - I28| for the period-two, synchronized states of a frac- 
tal growth network at level of construction L = 3, embedded 
in d = 2. (a) The upper curves correspond to the r.h.s of Eqs. 
I27ll - 128( l equal to —1 for both types of eigenvalues. Arrows 
indicate the boundary curve S3 (639) = +1. The interior re- 
gion bounded by these curves is where stable, synchronized, 
period-two states exist in the parameter plane (A, 7). (b) 
Magnification of the upper curves in (a), showing the gaps 
in the stability boundaries of the period-two, synchronized 
states. Arrows indicate boundary curves corresponding to 
several eigenvalues. The symbol (0) just beyond the bound- 
ary 5*3(636) = — 1 indicates the values of the parameters 7 
and A used in Fig. 6. 



VI. CONCLUSIONS 

The underlying inhomogeneous structure of fractal 
growth networks has significant effects on the spatial pat- 
terns that can be formed by reaction-diffusion processes 




FIG. 6: Inhomogeneous, period-4 state at parameter values 
7 = —0.296, A = 3.24. This pattern is a linear combination 
of the three eigenvectors associated to the eigenvalue 03 q . 



on these geometrical supports. In systems of interacting 
elements, such as the models considered in this article, 
the coupling matrix contains the connectivity of the net- 
work. The set of eigenvectors of the coupling matrix re- 
flect this connectivity. The spatial patterns that can arise 
in fractal growth networks are determined by the eigen- 
vectors of the diffusion coupling matrix M, which con- 
stitute a complete basis on which any spatial mode can 
be expressed as a linear combination. These eigenvec- 
tors have complex spatial forms but they are analogous 
to the Fourier eigenmodes arising in regular Euclidean 
lattices. On the other hand, the stability of the synchro- 
nized states is determined by the eigenvalues of M. The 
density distribution of these eigenvalues and their degen- 
eracy are nonuniform. These features affect the bifur- 
cation properties of dynamical systems such as coupled 
maps defined on fractal growth networks. The scaling 
structure of the synchronized, period-doubled states on 
the space of parameters of the system is similar for both 
uniform and fractal networks but the nature of the bi- 
furcation boundaries is different. For fractal growth net- 
works, the nonuniform distribution of eigenvalues leads 
to gaps in the boundary curves separating stable syn- 
chronized states that are not present for coupled maps 
on uniform lattices, where the spectrum of eigenvalues 
is continuous. These gaps allow for the selection of spe- 
cific spatial patterns arising from a uniform, synchronized 
state by appropriately changing the parameters of the 
system. 

Although we have presented only the simplest spa- 
tiotemporal patterns that can be formed on fractal 
growth structures, the formalism introduced in this ar- 
ticle can be applied to many other processes, such as 
excitable dynamics, nontrivial collective behavior, phase- 
ordering, domain segregation, turbulence, etc. The for- 
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malism is also useful for cellular automata models and 
continuous-time local dynamics on fractal growth net- 
works. 

The study of dynamical systems defined on nonuni- 
form spatial supports, such as fractal growth networks, 
provides insight into the relationship between topology 
and spatiotemporal phenomena in complex networks. 
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